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SUMMARY 


This paper presents the algorithms and logic used in the implementation 
of the microwave landing system back azimuth signals for the generation of 
navigation information for missed approach guidance. Navigation position 
estimates are based on range information from a randomly located DME and 
MLS back azimuth angular information. MLS volumetric coverage checks are 
performed to ensure that proper radio navigation inputs are being utilized. 
These algorithms and volumetric checks are designed so that they could be 
added to most existing area navigation systems (ARINC Characteristics 582 
and 583) with minimum software modification. 

INTRODUCTION 

The NASA TCV B-737 airplane has demonstrated the ability to navigate 
along a curved path approach to landing with an on-board navigation system 
utilizing position information from the time referenced scanning beam (TRSB) 
microwave landing system (MLS). While operating in the MLS environment, 
azimuth and vertical information was obtained from two scanning beams and 
distance information was provided by a precision MLS DME for use in the 
on-board navigation and autoland systems. The TRSB MLS provided azimuth 
information 60 degrees either side of the runway centerline, vertical 
information from the surface to 20 degrees elevation, and precision DME 
information 5 to 7 n.mi. from the MLS antenna. 

Missed approach navigation capability was recently added to the MLS 
system in the form of back azimuth and DME information. Hence, to fully 



utilize the TRSB MLS, it became necessary to incorporate the use of back 
azimuth information into the NASA TCV B-737 airplane navigation computer 
for use in generating missed approach guidance. 

The purpose of this report is to describe the equations and logic 
used to generate a navigation position estimate in the MLS back azimuth 
signal environment. The equations used to calculate position error components 
derived from 1) MLS back azimuth angular position measurements and randomly 
located DME range measurements and, 2) only MLS back azimuth angular 
position measurements are described. A summary of the NASA TCV B-737 
airplane's estimate process will also be described so that the total 
navigation process may be shown. 


SYMBOLS 

A - distance between the DME and MLS back azimuth antenna - m 
a - East component of A - in 

b - North component of A - m 

D - DME corrected for slant range - n.mi. 

D' - DME slant range - n.mi. 

DME - distance measuring equipment 

DF - position difference vector 

DP, - component of the position difference vector, DP, perpendicular 
-L to the runway centerline. 

DP^ - magnitude of vector DF^ - n.mi. 

dt - cycle interaction time of navigation computer - sec 

F - ellipticity constant, 0.003367 
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- altitude of aircraft above mean sea level - m 

- altitude of DME antenna above mean sea level - m 

- altitude of MLS back azimuth antenna above mean sea level - m 

- unit coordinate vectors 

- navigation update mode; inertial velocity, back azimuth & DME 

- navigation update mode; inertial velocity & back azimuth 

- position update gain 

- velocity update gain 

- outer radial limit of back azimuth volumetric check - n.mi. 

- inner radial limit of back azimuth volumetric check - m 

- axes of orthogonal coordinate system oriented towards True 

North 

- angle formed by the vector Z and a line between the airplane 

and DME antenna - deg m 

- Earth radius - m 

- Meridianal radius of curvature - m 

- Normal radius of curvature - m 

- unit vector perpendicular to the runway centerline 

- North & East components of system velocity estimate - KTS 

- North & East components of inertial ground speed - KTS 

- axes of orthogonal coordinate system oriented along the 

runway centerline 

- coordinates of vector Z transformed into the X', Y 1 coordinate 

system e 

- airplane's estimated position vector 

- magnitude of 1 Q the airplane's estimated position vector - n.mi. 

- North & East component of 1 the airplane's estimate position 

vector - n.mi. , 
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airplane's measured position vector 
magnitude of 2^ - n.mi. 

North & East components of 

vector of airplane's position estimated radially along thp 
measured azimuth angle 

magnitude of position vector Z R - n.mi. 

North & East components of position vector Z R 

angle formed by the DME antenna, origin, and measured 
airplane position - degrees 

North & East component of position estimate error - n.mi. 

North & East components of system velocity update - KTS 
latitude and longitude update estimate 
back azimuth angle relative to the runway centerline - deg 
relative angle between the DME & MLS back azimuth antenna - deg 
pi - 3.1416 

latitude and longitude of airplane position estimate - deg 
latitude, longitude of back azimuth antenna location - deg 
latitude, longitude of DME antenna location - deg 
runway heading to true North - deg 

vertical angular limit of back azimuth volumetric check - deg 
lateral angular limit of back azimuth volumetric check - deg 
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DISCUSSION 


NAVIGATION POSITION ESTIMATE DESCRIPTION 
Position Error 

The TCV B-737 navigation computer is software controlled to select and 
tune two appropriate DME and/or VOR stations in the vicinity of the airplane. 
The distance and/or azimuth information received from these stations is 
used to determine a position error from the previous position estimate. 

This position error is divided into North and East error components to be 
used in determining a new position estimate. 

Other sources of navigation information may also be used to determine 
position errors. While operating in the MLS back azimuth environment, 
these position error components are calculated with back azimuth and DME 
information. This estimate mode is shown on the pilot's electronic map 
display as IBD (inertial, back azimuth, DME). If no DME is available, the 
estimate mode is IBX. 


Position Estimate 

In the TCV B-737 airplane's navigation system, the position estimate 
calculations are used regardless of the manner in which the position error 
is determined. The first step in the position estimate process is to 
develop a system velocity update from the position error 

AV N = AV N + K 2 AP N 
A V E = AV e + K 2 AP e . 
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A system velocity estimate is obtained by summing the system velocity 


update with ground speed obtained from the inertial navigation system. 
v M = av m + V M 

N N N 


V = AV + V 
V E nv E V E 


A position update, in terms of latitude and longitude is obtained 
using the system velocity and the position error as follows: 

A<J> = (V N dt + K ] AP n )/R m 

A A. = (V £ dt + K ] AP e )/R n . 

This position update is based on an oblate spheroid earth model (ref. 1) 
by using appropriate radii of curvature in the North/South and East/West 
directions. 

R m = H a y G + R [;0 - 2F + 3F sin 2 4> e ) (North/South) 

2 

R n = H a / C + R E 0 + f sin 0) cos <f> e (East/West) 

The updated position estimate is found by summing the old previous 
position with the position update terms. 

<f> e = + A <f> 


A 

e 


= A + AA 

e 


BACK AZIMUTH/DME (IBD) POSITION ESTIMATE 
General Solution 

Figure 1 shows the geometry of the MLS back azimuth antenna, the DME 
antenna, the airplane's position estimate and its position as measured with 
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back azimuth and DME information, and the runway with an extended centerline. 
An orthogonal coordinate system, with its origin placed on the back azimuth 
antenna, is oriented in a True North direction. Since the DME may be either 
randomly or manually selected, it is not required to be co-located with the 
back azimuth antenna. 

By determining the relative geometry between the back azimuth antenna 
and the DME and the airplane's position measured with the back azimuth 
and DME information, the vector between the origin and the airplane's 
measured position, Z^, may be found. The vector between the origin and the 
airplane's previously estimated position, Z g , may be found directly from 
their known latitudes and longitudes. The desired position error used to 
form a new position estimate may then be found by subtracting from Z^. 

Calculation of the Airplane Position Vector 

Figure 1 shows that is one side of the triangle formed between the 
origin (back azimuth antenna), the airplane's measured position, and the 
DME antenna. Known quantities used to determine Z M in this triangle include 
the runway heading, ij> R ; the slant range of the DME, D ' ; the back azimuth 
angle relative to the runway centerline, n; and the latitudes and longitudes 
of the DME antenna, <j>^£, anc * the back azimuth antenna location 



Figure 2 shows the angular geometry and distance between the origin 
and DME antenna, A. The length A and the relative angle between the DME 
antenna from the North axis, y, remains constant and is calculated only once. 
If a new DME is tuned, then A and y are recalculated. 
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The length of A is determined by vectorial ly summing its components 


a and b. 


^ X DME " V< 6 °) cos 
<+DME - *o )(60) 


♦o + ^DME 


2 2 
a^ + o 


The angle y is found by 
.. _ *,--1 fa 


y = tan 


N 


The angle a, formed by side A and vector Z^, may continuously vary and 
is calculated 20 times per second. Angle a is found by taking the absolute 
value of the sum of the azimuth deviation angle, and the angle y, and the 
difference of the runway heading. 


a = y 


- ton - n) 


If the absolute value of a is greater than n, then: 
a = 2JI - a 

The magnitude of the DME reading as measured in the airplane, D 1 , is 
the slant range distance between the ground-based DME antenna and the 
airplane. This distance was slant-range-corrected to determine the ground 
distance, D, between the airplane and DME 


D = D sin cos' 


Vc ~ DM E 
D* 


Angle P is formed by the vector and side D. This angle may 
continuously vary and is calculated 20 times per second. Knowing the 


relation 


P = D , 

sinP sina 



angle P is determined by 

n • -1 r A . 

P = sin p- sina 

It must be determined if angle P is an oblique or an acute angle. This 

is accomplished by comparing the square of side A with the sum of the square 

of side D and the square of the magnitude of the estimated airplane position 

vector, T . If 

e P 

2 ? c 

A > D + Z g , then P is oblique and set equal to 
P = tt - P. Otherwise, P = P. 

The magnitude of Z £ is used as an approximation for 2^. The magnitude 

of Z is: 
e 

Z = 
e 

in which 




(4> e - 4> Q )(60) 

( A e ' V (60) cos 


*o + * 


e 


2 


The magnitude of Z^ is found by 

Z^ = A cos (a) + D cos (P). 

The North and East components of Z^ are found knowing the angle, 
(i|> R -n)> between the North axis and Z^. 

Z MN = Z M C0S ^R" n ) 

Z ME = Z M Sin ^R' n ) 


Calculation of the Desired Position Error in North & East Components: AP^, 

AP^-IBD Update Mode 

The desired position error, in North and East components, is now found 
by subtracting the North component of 1 Q from the North component of and 
the East component of Z g from the East component of Z^, 
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AP N Z eN ' Z MN 


AP C = Z c 
E eE 


-ME 


These position error components are then used directly in the 
navigation position estimate algorithms. 


BACK AZIMUTH ONLY (IBX) POSITION ESTIMATE 
General Solution 

In the event that a DME can not be tuned, the navigation computer will 
utilize the back azimuth signal and inertial velocity to determine a position 
estimate. Figure 3 shows the geometry of the back azimuth antenna, the 
airplane's position estimate and its estimated position on a measured back 
azimuth angle, and the runway with an extended centerline. An orthogonal 
coordinate system, with its origin placed on the back azimuth antenna, is 
oriented in a True North direction. 

Since no DME information is available, no radio position errors can be 
developed in a radial direction from the back azimuth antenna. However, 
inertial velocity is still utilized in the radial direction and will supply 
inputs for a new position estimate in the position estimate algorithms. 

Radio position error in the IBX update mode is limited to a direction 
perpendicular to the runway centerline. 

North and East position estimates are found in the following manner. 

A position difference vector, DP, is found by subtracting the estimated 
position vector, Z g , from the estimated position vector on the measured 
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back azimuth angle, Z R . The component of DP perpendicular to the runway is 
found and broken into North and East components. These components are used 
in the navigation position estimate algorithms. This process is repeated 
20 times per second. 

Calculation of DP 

To determine DP it is necessary to calculate the position estimate 
vector, Z g , in North and East components from the latitudes and longitudes 
of the last position estimate and the back azimuth antenna location. 


'eN 


(60) (0> e - # 0 ) 


^ + <t>. 


Z eE = (60) (A e - A q ) cos 

/S. /N 

7 = Z i + Z j 

e eE eN ° 

A vector Z R , of the airplane's estimated position along a measured 
back azimuth angle is determined in the following manner. Since no radio 
updates can be obtained in a radial direction from the back azimuth antenna, 
it will be assumed that the estimated radial distance from the origin is 
correct. Hence, the vector length of Z D and Z will be the same (directions 
may differ to obtain DP). 

The length of Z D and is: 

K £ 


Z R ■ Z e =1 Z e N " ♦ z eE Z 

The North and East components of Z R are found knowing the angle, 
(i|j R -ri), between the vector Z R and the North axis. 


Z RN 

Z R 

cos 

(^-n) 

Z RE = 

Z R 

sin 

(^ R -n) 

Z 

z 

i + 

Z„., j 

R 

RE 


RN J 


11 



DP is found vectorial ly 



Calculation of the Desired Position Error: AP^, AP £ - IBX Update mode 

The magnitude of the component of DP perpendicular to the runway 
centerline is obtained by the vector dot product of DF with a unit vector, 
u, perpendicular to the runway centerline. This results in the magnitude 
of the desired position error. This magnitude is then multiplied times the 
unit vector to obtain the desired position error in North and East components. 
These components are then used directly in the navigation position estimate 
algorithms. 

The unit vector, u, shown in figure 3, is: 

A A A 

u = - cos (ip R ) i + sin (ip R ) j 

The magnitude of the component of DP perpendicular to the runway 
centerline is: 

DP, = DP ' u 
1 

DP 1 ' ' (Z RE - Z eE> C0S <*R> + (Z RN ' Z e N> s1n 'V 

The North and East components of the desired position error are: 

AP n = DP^ sin (ip R ) 

AP £ = - DP^ cos (ip R ) . 

These position error components are then used directly in the 
navigation position estimate algorithms. 


12 



BACK AZIMUTH/AIRPLANE POSITION VALIDITY CHECK 

Since random DME selection, automatic frequency tunning, and other means 
of automatic software control will be utilized in the navigation computer, 
a check must be made to ensure that the navigation computer is using the 
appropriate navigation data for the area in which the airplane is flying. 
Obviously, if improper navigation data is being utilized, position estimates 
based on that data must be inhibited. 

To preclude the possibility of using the wrong MLS back azimuth 
navigation data, a volumetric geometry check is made to determine if the 
airplane's position estimate is within the MLS back azimuth boundaries. 

If the airplane is not within these boundaries then MLS updating is 
inhibited. These boundaries include lateral, radial, and vertical limits 
of coverage as shown in figure 4. The lateral limit of the back azimuth 
angle of coverage is ±u> from the back azimuth antenna. Radial limits 
require that the airplane be within L n.mi. of the back azimuth antenna, 
but not closer than M meters. 

A new orthogonal coordinate system, with its origin located at the 
back azimuth antenna and x 1 axis parallel to the runway centerline, was 
used to make the geometric checks. To make the checks, the North 
and East components of the estimated position vector Z , are transformed 
into the new x', y' coordinate system as follows: 
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X 'e = Z eN cos (*r> + Z eE s1n 

V 'e = Z eN sin (*r) - Z eE cos <%>• 

The lateral azimuth check is then 

-X' e tan ui < Y' e < X' g tan w. 

The radial check is 

M meters < X' < L n.mi. 
e - 

The vertical check is 

H a/c - H 0 + X ' e tan 

CONCLUDING REMARKS 

Area navigation systems may use various radio inputs to calculate a 
position estimate for course guidance. Flight tests of the NASA TCV B-737 
airplane using the Time Referenced Scanning Beam microwave landing system 
have shown that an area navigation system may utilize front course signals 
to provide guidance for curved path approaches to landing. 

MLS back azimuth signal coverage has been added to the microwave 
landing system to provide guidance during missed approaches. These back 
azimuth signals are utilized on the TCV B-737 airplane by algorithms that 
generate position error information to be used in the area navigations 
system. These algorithms are designed so that they may be added to an 
existing area navigation system (ARINC Characteristic 582 and the more 
sophisticated systems under ARINC Characteristics 583) with minimun 
software modification. 
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Figure 1.- Back azimuth, DME, airplane estimated and measured positions, and runway with an 

extended centerline geometry - IBD update mode. 
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